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ABSTRACT 

Globular cluster systems of 17 elliptical galaxies have been studied in the 
Coma cluster of galaxies. Surface-brightness fluctuations have been used to de¬ 
termine total populations of globular clusters and specific frequency (Sn) has 
been evaluated for each individual galaxy. Enormous differences in Sn between 
similar galaxies are found. In particular, Sn results vary by an order of magni¬ 
tude from galaxy to galaxy. Extreme cases are the following: a) at the lower end 
of the range, NGC 4673 has Sn = 1.0 ± 0.4, a surprising value for an elliptical 
galaxy, but typical for spiral and irregular galaxies; b) at the upper extreme, 
MCG +5 -31 -063 has S N = 13.0 ± 4.2 and IC 4051 S N = 12.7 ± 3.2, and 
are more likely to belong to supergiant cD galaxies than to “normal” elliptical 
galaxies. Furthermore, NGC 4874, the central supergiant cD galaxy of the Coma 
cluster, also exhibits a relatively high specific frequency (Sn = 9.0 ± 2.2). The 
other galaxies studied have Sn in the range [2, 7], the mean value being Sn = 5.1. 
No single scenario seems to account for the observed specific frequencies, so the 
history of each galaxy must be deduced individually by suitably combining the 
different models (in situ, mergers, and accretions). The possibility that Coma 
is formed by several subgroups is also considered. If only the galaxies of the 
main subgroup defined by Gurzadyan & Mazure (2001) are used, a trend in Sn 
arises in the sense of Sn being bigger in higher density regions. This result needs 
further confirmation. 

Subject headings: galaxies:clusters:individual (Coma)-galaxies:star clusters 
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1. Introduction 

Globular clusters (GCs) are thought to be among the oldest objects in the Universe; so, 
although they represent only a small fraction of the total luminosity of the galaxy hosting 
them, they provide useful information about the galaxy formation process and its initial 
evolutionary stages. By studying globular cluster systems (GCS), the memory of the host 
galaxy is probed. 

A fundamental, widely used parameter describing a GCS is the specific frequency, Sn, 
introduced by Harris & van den Berg (1981). Sn is the number of GCs normalized to an 
absolute magnitude for the parent galaxy of My OT = —15 mag: 

S N = Nqq io a4(M v OT + 15) , (1) 

where N tot is the total number of GCs in the galaxy. Sn gives the number of GCs per unit 
luminosity of the host galaxy. Spiral and irregular galaxies have typical values of Sn < 1 
(van der Berg & Harris 1982). Dwarf (Durrell et al. 1996; Miller et al. 1998), elliptical 
(Harris & van den Berg 1981), and lenticular galaxies (Kundu & Whitmore 1998; Chapelon 
et al. 1999) have Sn between 2 and 6, where larger values seem to correspond to galaxies 
located in high-density environments (Harris 1991). Giant cDs located near the centers of 
rich clusters of galaxies have the largest specific frequency, Sn ~ 10 — 20 (Harris & van den 
Berg 1981; Harris, Pritchet, & McClure 1995; Blakeslee & Tonry 1995; Bridges et al. 1996), 
with some exceptions (Kaisler et al. 1996). 

The origin of the observed values of Sn is not clear. A basic question is whether 
Sn depends on the properties of the host galaxy and/or the environment. If only elliptical 
galaxies are considered, environmental density seems to play a fundamental role in Sn (Harris 
1991). These results, which are based on the compilation of several studies of elliptical 
galaxies closer than the Virgo cluster, suggest that the formation of GCs has been two or 
three times more efficient in rich environments than in poor ones. On the other hand, the 
luminosity of the host galaxy also has an influence on the properties of the GCS (Harris, 
Harris, & McLaughlin 1998). 

Forte, Martinez, & Muzzio (1982) and Muzzio (1987) have suggested that in clusters 
of galaxies a population of intergalactic GCs (IGCs) should exist. A galaxy could enlarge 
its GC population via gravitational capture of IGCs. If they are distributed according to 
the gravitational potential of the whole galaxy cluster, then the influence of the environment 
on Sn is justified (White 1987; West et al. 1995). On the other hand, Blakeslee (1997) 
suggested that Sn could be strongly influenced by gas removal by tidal disruption during the 
formation of the cluster. This idea was developed further by Blakeslee, Tonry, & Metzger 
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(1997). The star formation would be inhibited as a result of the gas loss, increasing Sn- 
Moreover, Harris et al. (1998) argued that explaining the observed Sn in M87 via gravita¬ 
tional capture of IGCs, requires assuming an unrealistic spatial distribution for intergalactic 
GCs. 

The origin of the GCSs themselves in giant elliptical galaxies is still not understood. One 
of the main problems in understanding this phenomenon is the lack of data. The difficulty 
in observing distant GCSs is due to the low resolution of available telescopes and to the 
low brightness of the individual GCs, which are too faint to be detected even at moderate 
distances. At ~ 100 Mpc (the distance of the Coma galaxy cluster), only the brightest GCs 
are visible from ground-based observations (Harris (1987) and Thompson & Valdes (1987) 
performed the first studies of GCSs in Coma); high-resolution telescopes (such as the Hubble 
Space Telescope ( HST ), see for example Harris et al. (2000), Kavelaars et al. (2000)) and long 
exposure times are required to detect them directly. But the surface-brightness fluctuations 
(SBF) technique, which we use in this paper, allows us to detect GCSs in elliptical galaxies at 
distances of the order of 100 Mpc, using small to intermediate-size, ground-based telescopes 
and reasonable exposure times (Blakeslee & Tonry 1995). 

The SBF technique was introduced by Tonry & Schneider (1988) with the aim of 
measuring distances. Comparing the fluctuation signals produced by the stellar population 
of a galaxy with SBF measurements of nearby galaxies which have externally calibrated 
distances, accurate estimates of distances can be obtained up to ~ 40 Mpc (Tonry et al. 
2000, 2001). For smaller distances, undetected GCs contribute as a perturbation of the total 
fluctuation signal of the galaxy. But in the case of more distant galaxies the fluctuation 
signal produced by the GCs, which are mostly unresolved, dominates and the contribution 
of stars becomes negligible (Wing et al. 1995; Blakeslee & Tonry 1995; Blakeslee, Tonry, & 
Metzger 1997; Blakeslee 1999). 

In this paper, we use the SBF technique to evaluate the GCS properties of elliptical 
galaxies in the Coma cluster. Coma is a good laboratory for studying GCSs, because it is a 
rich cluster and its ellipticals have a wide range of luminosities and environmental densities. 
After describing the SBF technique and its use in obtaining GCS properties (§3), we present 
a test for evaluating the capacity of the technique (§4), which we then apply to the Coma 
cluster galaxies (§5). 
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2. Observations and Data Reduction 

This study is based on observations done in 2000 April 25 and 27 with the 2.5 m INT 
at the Roque de los Muchachos Observatory on La Palma, using the Wide Field Camera, 
which consists of four 2048 x 4100 pixel chips with a scale of 0.333 "/pix, providing a held 
of view per chip of 11.38' x 22.75'. The total held of view is ~ 34' on a side. 

In Table 1, the location of the target helds and all the galaxies studied in each held 
are listed. Column 1 gives the name adopted in this article for each target held, columns 
2 and 3 give the coordinates of the held and column 4 gives the galaxies analyzed in each 
held. Twelve exposures of 300 s with the Sloan R Filter were taken for each held, except for 
the central held of Coma (RC-1) where thirteen exposures of 300 s were taken; this resulted 
in a total integration time of 3600 s (or 3900 s for RC-1). In Table 2, a summary of the 
observations is given. Column 1 lists the held names, column 2 the date of the observation, 
columns 3 and 4 give the exposure times and the used hlter, and column 5 gives the seeing, 
which is between 1" and 1.2" in all cases. The bias level was subtracted from the individual 
frames using the overscan region of the CCD. The data were then battened with twilight 
sky hats. The hnal image was produced by calibrating the individual frames, dithering 
between exposures and calculating an average of the frames. In Fig. 1, the central chip of 
the CCD is shown for the target held RC-1, where nine studied galaxies are marked. This 
held corresponds to the center of the Coma cluster. 

Photometric calibration was performed for the second night. The observations of the hrst 
night were then calibrated comparing the held RC-1 observed in both nights. The extinction 
and photometric calibration were calculated using hve measurements of 23 Landolt (1992) 
standard stars. The equation transforming instrumental to standard magnitudes is 

R — r + AX = a,i + bi x (V — R), (2) 

where R is the magnitude in the standard system, r is the instrumental magnitude, X is 
the airmass (all the observations of target objects have airmass smaller than 1.1) and A is 
the extinction, for which a value A = 0.145 was obtained. (V — R) is the colour index of 
the program objects. Since no observations were done in the V Filter, a typical value for 
the colour index of GCs of (V — R) — 0.46 was adopted (Ajhar et al. 1994). a* and A are 
the zero point and the color term for each chip (i=l, 2, 3 and 4). Our results for a, and 
bi are listed in Table 3. The sigma of the photometric calibration varies between 0.01 and 
0.02 mag. Calibration results for each individual galaxy are shown in Table 4, where rri* is 
the magnitude of an object yielding one unit of hux per unit time. This magnitude will be 
used later. The sigma of rn\ represents the sigma of the photometric calibration. Finally, as 
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Coma is nearly perpendicular to the Galactic plane, interstellar extinction is neglected. 


3. Surface Brightness Fluctuations: the Technique 

The concept of surface-brightness fluctuations was introduced by Tonry & Schneider 
(1988), who noted that in the surface photometry of a galaxy far enough to remain unre¬ 
solved, a pixcl-to-pixel fluctuation is observed due to the Poisson statistics of the spatial 
distribution of stars, GCs, background galaxies, foreground stars, etc. This is correct for 
a galaxy with no absorption or emission due to dust or gas; so this technique would only 
be applicable to early-type galaxies. The variance of the fluctuation depends on the stellar 
population, the GCS, background galaxies, foreground star, and, of course, the distance. 


3.1. Sources of Variance: GCS Properties 

The following treatment has been described by Blakeslee & Tonry (1995). It can be 
assumed that the total pixcl-to-pixel variance of an image (<j^ OT ) is the sum of all the 
independent contributions due to the different components in the image, such as the stellar 
population (crg P ), GCs (ctq C ), background galaxies (<v| G ), foreground stars (crf s ), read-out 
noise (cry) and photon shot noise (<7p h ): 

2 2 i 2 i 2 i 2 i 2 i 2 /o\ 

^TOT = Gp + GjC + ff BG + G S + Go + %h ( 3 ) 

The pixel-to-pixel variance produced by a class of object can be evaluated as the second 
moment of the luminosity function of that object population. In general, the brightest 
individuals could be resolved in the image. For the analysis of the SBF, these objects must 
be masked. If all sources brighter than a limiting flux (./’u m ) are masked in the image, then 
the variance produced by the remaining part of the population is: 


°pop — 


n(f)f 2 df, 


( 4 ) 


Jo 

where n(f ) is the luminosity function; i.e., the number of population objects per unit flux 
per pixel. We can put this equation in terms of magnitudes via the relationship 


m = —2.5 log (/) + ml, 


( 5 ) 
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/ being the flux (e s 1 pix and m\ the magnitude of an object yielding one unit of flux 
per unit time. 

If n(f) is known, the variance produced by the non-masked part of the population can 
be evaluated. In the case of GCs, a Gaussian shape for the luminosity function (GCLF) can 
be assumed (Harris 1991): 


( \ b(C —— mP)^ / 2<r^ 

n GC {m) = —f=^e 1 > ' , 

v27ra 


( 6 ) 


where Nqc is the total number of GCs per pixel. Substituting eq. 6 into eq. 4, the variance 
due to the non-masked globular clusters can be expressed as: 


a 2 c = J^GC_ 1 Q 0.8m* f e -(m-mO)V2^ e -0.81n(10)m dmj (7) 

V 2n <j J m c 

where m c is the magnitude corresponding to the limiting flux /i; m . 

It can be seen that the variance produced by the GCs is proportional to their total 
surface density. In this way, assuming values for m° and cr in the GCLF, and determining 
the variance produced by GCs, <Jq C from the image analysis, their surface density can be 
obtained. In galaxies with large angular size, the spatial structure of the GCS can be also 
evaluated by dividing the galaxy into rings and measuring the GCs’ variance in each one. 

What is directly obtained from the analysis of the image is ct pot . Hence a treatment 
of the remaining sources of variance is required. First, at the distance of Coma, using eq. 7 
and with the help of the Padua stellar evolution library (see Bertclli et al. (1994)), it can 
be shown that the ratio Ug P /<jQ C is nearly zero, so <jg P can be neglected in eq. 3. Also, as 
Coma is nearly perpendicular to the Galactic plane, foreground stars are very few, and cr| s 
can be neglected. Hence, 


°TOT ~ ^GC + °BG + a ph + °fo• (8) 

As we will see later, the SBF technique involves spectral analysis of the signal. This 
provides two quantities: P i and P o- P o is the total PSF-convolved variance, produced by 
all objects whose spatial flux distribution is convolved with the PSF. P\ is the non-PSF- 
convolved variance. As a result, 


P 0 ~ (Tg C + UgQ 


(9) 
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and 

P 1 = ^ro + ^ph- ( 10 ) 

Consequently, the only contribution to we have to be concerned with is er| G . Once this 
is removed, <t gc will be obtained from Pq. 

The following law can be assumed for the luminosity function of background galaxies 
(BGLF): 


n BG (m) =Ax N bg l(P m , (11) 

where A is a constant for normalization, TVbg is the total number of galaxies per pixel and 
7 is the slope of the magnitude distribution. Tyson (1988) found 7 ~ 0.39 in the R filter. 
Using eq. 4, we find 


a 2 BG = A * Nbg 10 7mc 10 o.8K- mc ) 

BG In 10 ( 0 . 8 - 7 ) 

Once the parameters A, N BG and m c are estimated from an analysis of the image, as 
we will see later, and rn\ from the photometric calibration, cr| G can be calculated. 



3.2. Measuring er GC 

As we have seen, because stars, background galaxies, and GCs are convolved with the 
PSF, while the read-out and photon shot noise are not, if a power spectrum of the image 
is computed, the non-convolved variance can be separated from the remaining terms. The 
power spectrum of the image, P{k), has the form 


P(k) = P l + P 0 E(k), (13) 

where E(k) is the power spectrum of the PSF convolved with the window function (which 
is unity in all the regions we want to study and zero in the rest of the image), P 0 is the 
PSF-convolved variance (P 0 = Ugp+a GC + < 7 BG+ cr FS+'")> an( l Pi is th e non-PSF-convolved 
variance (Pi = ofo+°"ph)- 

Before performing the power spectrum of the image, a preliminary image treatment has 
to be done. 



3.2.1. Preliminary image treatment 


After the initial image reduction was done, and because there are some galaxies located 
in each chip of the CCD mosaic, different regions of the chips were selected, each containing 
a galaxy (see Fig. 1) that was treated individually. In all cases, a smooth fit of the galaxy 
profile was performed using the IRAF package ISOPHOTE. An iterative method was used: 
first, bright objects that could affect the galaxy profile fit were masked. The task ELLIPSE 
was used to fit this profile, and the result was subtracted from the original image in order 
to obtain a map of the bright objects alone. These were then subtracted from the initial 
image and a second fit of the galaxy profile was made. This procedure was repeated until 
the precision of the final fit was good enough. 

Once the final galaxy profile was performed and subtracted, a residual zero mean image 
is obtained where the PSF-convolved variance can be measured. In this image the power 
spectrum is only perturbed at low wave numbers (Jensen, Tonry, & Luppino 1999), so this 
part of the power spectrum will not be taken into account in our analysis. 

Using the DAOPHOT task DAOFIND, all the objects 4<r above the zero mean of the 
residual image are selected. Those brighter than some limiting magnitude (m c ) were masked 
out creating a window function whose value is zero where the bright object is located and 
unity in the rest of the image. The magnitude m c must be chosen appropriately. Note that 
if a very bright m c is used, Po is dominated by the non-masked brightest objects (since it 
is the second moment of the luminosity function of the population which is at work). In 
our study, brightest objects are background galaxies, so the GC contribution is masked out. 
Because the CG contribution will be evaluated, as many bright objects as possible must 
be masked out. On the other hand, if m c is chosen too faint, the results could be affected 
by incompleteness, resulting in an incorrect estimate of the unresolved background galaxy 
contribution to Pq. An analysis of this question has to be made in each case. In this article, 
we choose m c = 23.5. 

In large-angular-size galaxies, the spatial distribution of the GCS can be analyzed. In 
those cases, the set of ring-shaped regions in which the Pq will be measured can be defined 
making use of the galaxy profile fit. We can select different isophotes to define the annular 
sections, in this way, elliptical annular regions are obtained at different distances from the 
center of the galaxy. Multiplying the previous window function with each defined annular 
region, a set of new window functions is obtained, one for each ring that we shall be studying. 
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3.2.2. Measuring <7 q C 

Multiplying the window function of a selected region by the residual image, the residual 
masked image of the region is obtained. After this, the power spectrum of the residual 
masked image is computed. The result is a two-dimensional power spectrum, which is 
radially averaged in order to obtain the one-dimensional power spectrum of the selected 
region. This has the form expressed in equation 13. In order to evaluate Pq, the power 
spectrum of the PSF convolved with the window function of the region is needed. Bnt E(k) 
can be approximated by the power spectrum of the PSF alone, Ppsf(^)> with a negligible 
error (Jensen et al. 1998), so eq. 13 transforms into 


P(k) = P 1 + P 0 P PSF (k). (14) 

Using an empirical form for the PSF, obtained from some bright, isolated stars in the 
same image being analyzed, and fitting the power spectrum of the images with eq. 14, the 
quantities Pq and Pi can be evaluated. In large-angular-size galaxies, the spatial structure of 
the GCS can be deduced from the radial dependence of P 0 . In the case of smaller galaxies, 
where we cannot measure radial dependences, only the total number of GCs is evaluated. 


4. Testing the Technique Using Synthetic Images 

Putting into practice what we have presented in the preceding sections, in this section 
we will make a consistency test to check our procedure for measuring the power spectrum 
normalization. The influence of the noise level in the images is also tested. Synthetic images 
were constructed with this purpose in the following manner. The magnitudes of the GC were 
randomly calculated according to a Gaussian distribution. The parameters of the GCLF were 
a = 1.4 and m° = 20. This is equivalent to supposing that the GCS is at a distance of 4 Mpc. 
At this distance, GCs have a stellar appearance, so simulating them with an empirical PSF 
is adequate. Using an empirical PSF and the DAOPHOT task ADDSTAR, we generated 
10 000 artificial GCs uniformly distributed across a 1000 2 pixels image. Six different values 
for the non-convolved variance were used (Fig. 2) and the SBF technique was applied to the 
final synthetic images in order to recover the total number of GCs. Since all the objects in 
the image are GCs, it is not necessary to mask the brightest sources, so we used m c = — oo 
and m* = 25 in equation 7. In this way we obtain ctq C = 0.278./Vtot, where Ntot is the 
estimate for total number of GCs in the synthetic image. Finally, the SBF technique will 
produce for each image the values of Pi=cr 2 oise and P 0 =°gc- 
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In Fig. 3, the power spectrum of the synthetic images for the six different noise levels 
are presented. Table 5 lists the properties of the synthetic images (columns 1 to 3) and the 
results of the SBF analysis (columns 4 to 6 ). The test goes from images with no noise at 
all (image a) to completely noise-dominated images (images e and /). Here, a 2 oise is the 
input noise level, Nq C is the number of GCs introduced in the image, Pq and P\ are the SBF 
results for each image, and Nqq is the total population of GCs estimated using equation 7. 

It is very interesting, even in the images with extremely high noise variance, that SBF 
produces good results. In all cases, there is excellent agreement between the parameters of 
the input images and those recovered by the SBF technique, both for <r 2 oise and ctq C . We 
only observe a small systematic excess in the estimated number of GCs. This excess could 
be due to small differences between the real input GCLF (which includes random effects) 
and the exact Gaussian-shaped GCLF used to evaluate the final result in equation 7. In 
order to check this, we calculated < 7 q C directly from the magnitudes of all the 10 000 GCs 
added to the images in the following way: 

10 000 10 000 

4c =£/< = £ (15) 

i=l i= 1 

The result was (Tq C = 2926 (e _ /pixel) 2 , in full agreement with the obtained SBF results 
(column 5 in Table 5). With the former test we have shown that SBF is a very powerful, 
selfconsistent technique, and that very good results can be obtained even with completely 
noise-dominated images, where traditional techniques do not detect any object. 


5. Results 

5.1. Measuring the BGLF and <t 2 g 

The following treatment is the same for all studied galaxies, but with the aim of illus¬ 
trating the methodology, we first present the procedure for the galaxy NGC 4874 in detail 
and then the results for the rest of galaxies. 

Once the residual image of NGC 4874 is created, cosmic-ray events and bad pixels are 
masked. In order to measure < 7 b G , we must detect all objects brighter than 4 a above the 
zero mean SBF. Using the DAOPHOT task DAOFIND, 649 objects were detected. The 
photometry was made using ALLSTAR. In Fig. 4, the luminosity function of all detected 
objects is shown. The solid line represents the BGLF (eq. 11) with 7 = 0.39 (Tyson 1988) 
and scaled to our counts. The BGLF in the case of NGC 4874 is: 
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n Ba (m) = 3.17 x lCr 13 10«»" ( galaJd<B ) . (16) 

\ mag x pix J 

The fit is very good, so we can conclude that almost all point sources are background 
galaxies, as expected. 

Now, <r| G can be estimated using this BGLF and eq. 12. The uncertainty in the 
background variance must be evaluated taking into account the uncertainty of the BGLF 
slope as well as possible departures from a pure power-law model. Recent determinations 
of the i?-band faint galaxy luminosity function result in slopes in the range 7 = 0.30 — 0.34 
(Steidel & Hamilton 1993). Considering the result by Tyson (1988), 7 = 0.39, we have 
assumed that 7 = 0.30 — 0.40 is a reasonable range for the slope value. We have done 
simulations using these two extreme values, resulting in a 20 % variation in the background 
variance. We assume that this 20% is the uncertainty in the background variance. In 
particular, in the case of NGC 4874, for which m\ = 24.44 (see Table 4) and m c = 23.5, we 
obtain: 


^Ig — (29.4 ± 5.9)1(T 4 (—) . (17) 

\sx pix / 

This procedure is repeated for each galaxy. In Fig. 5, the luminosity functions of all 
detected objects are presented for each galaxy. Again, the solid line represent the BGLF 
with 7 = 0.39, scaled to our counts. Using these BGLFs, (j| G are obtained for each galaxy. 
Results are listed in Table 6 . 


5.2. Measuring errand the GCS properties 

As in §5.1, we first describe in detail the procedure used to measure <r GG for NGC 4874 
and give the results for the remaining galaxies. 

Before computing the power spectrum, all objects brighter than m c = 23.5 are masked 
and the ring-shaped regions are created in the galaxy for the analysis of the spatial distri¬ 
bution of GCs. The characteristics of all the regions are given in Table 7 (columns 1 to 3). 
t min and r max are the minimum and maximum distances of the selected ring-shaped region 
to the galaxy center, and A is the area of this region. 

The power spectrum of all the masked regions is then obtained and the SBF evaluated. 
In Fig. 6 , the SBF fit is shown for all the eight ring-shaped regions studied in NGC 4874. 
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Spatial variations of the PSF along the CCD mosaic may introduce a significant uncer¬ 
tainty in Pq. To limit this efect, we have used a PSF template for each of the four chips 
of the WFC. Moreover, in order to evaluate the uncertainty produced along a single chip, 
SBF fit was repeated for NGC 4874 using eleven PSF templates, which were computed using 
eleven PSF stars located at different positions of the chip. The diferences in Pq from fit to fit 
were less than 2%. This 2% was added in quadrature to the sigma of the P 0 measurements 
of all galaxies. The results for P 0 are listed in column 4 of Table 7. 

Subtracting < 7 g G from P 0 in all regions, cr GC is obtained. The results for a GC are also 
given in Table 7 (column 5). Note that NGC 4886 is located in the NGC 4889 ring-shaped 
region with [r mm , r max ] = [71.9", 104.7"], so a GC of NGC 4889 in this region has to be also 
subtracted to the NGC 4886 P 0 results. 

The GCLF parameters must be introduced in eq. 7 to obtain the number of GCs from 
Oq C . The GCLF for giant ellipticals seems to be universal. There is no empirical evidence 
showing galaxy-to-galaxy differences in GCLF parameters (Harris et al. 1998). Moreover, 
Kundu et al. (1999) studied 1057 GCs in the inner region of M87 using HST observations. 
This work shows no evidence for radial variations in the GCLF. Kavelaars et al. (2000) found 
rriy = 27.88 ± 0.12 for NGC 4874. This value comes from a constrained fit with cr = 1.40. 
Assuming a mean value (R — I) = 0.55 and (V — I) = 1.01 for GCs (Ajhar, Blakeslee, & 
Tonry 1994), m° R = 27.42 ± 0.12. 

Kavelaars et al. (2000), based on the reviews of Whitmore (1996) and Harris (1999), 
quote estimates of the intrinsic galaxy-to-galaxy scatters in a and m° R of 0.05 and 0.15 mag, 
respectively. Althought, from Ferrarese et al. (2000), one might estimate uncertainties twice 
as large. 

In order to include the influence of galaxy-to-galaxy variations in both a and m° R in 
our results, we adopted a = 1.40 ± 0.05 and m R = 27.42 ± 0.20 as our fiducial values and 
tested the influence that 1-sigma variations in both a and m° R have on <j gc /Nqc- To this 
purpose, eq. 7 was used to evaluate a GC using as input parameters the following sets of (a, 
m° R ): (27.22, 1.40), (27.62, 1.40), (24.42, 1.35) and (27.42, 1.45). The adopted results for 
(Xqc/Agc will be the mean ± a of the sets’ results. 

For NGC 4874 the next relation is obtained: 


a 2 GC = (0.067 ±0.016) x N GC 


s x pix 


(18) 


This calibration depends on the photometric calibration of each target held, so it is 
different for each galaxy. All a GC /Noc calibrations are presented in Table 8, where the ratio 



-13- 


°gc/ n gc is given for the different sets of imput parameters (columns 2-5) and for the final 
adopted value (column 6). 

Coming back to NGC 4874, eq. 18 can be used to obtain Nqc in each region, and 
multiplying by the area of the selected region, the number Vqq 1011 of GCs can be derived (see 
Table 7, column 7). If we consider the whole galaxy, the total population of GCs is obtained. 
For NGC 4874, we obtain 17600 ±4100 GCs. This is the number of GCs in the whole galaxy 
except the central region, which could not be analyzed due to the difficulty of making a good 
galaxy profile fit near the center. If the total population of GCs is sought, an estimate of the 
number of GCs in the central region is required. In this work, we assume that the surface 
density of GCs in the central region is the same as that obtained in the inner ring-shaped 
region analyzed in the galaxy. This approximation produces an underestimate of the real 
number of GCs in the center of the galaxy, but as the size of the region not studied is very 
small (in some cases negligible), the error introduced is also small. In most galaxies, the 
estimated number of GCs in the central region is about 5-10% of the total population. 

Extrapolating to the central region of the galaxy, we obtain a total population of 18200± 
4100 GCs for NGC 4874. 

In Fig. 7, the radial dependence of the GC surface density is presented. Open dots 
represent our results obtained with the SBF technique, and filled dots represent the results 
obtained using the HST and traditional techniques (Harris et al. 2000). Interestingly, when 
analyzed with the SBF technique, ground-based data produce similar results to those of the 
HST except for the galactic halo, where the surface density of GCs obtained in this work 
is bigger than that obtained by Harris et al. (2000). In particular, the total population 
obtained by Harris et al. (2000) is 9200 ± 1500 GCs. The result of the present study is closer 
to that in Blakeslee & Tonry (1995), who obtained a total number of GCs of 17260 ± 2030. 

The total number of GCs can be used to derive the specific frequency, 5%. The RC3 
(de Vaucouleurs et al. 1991) catalogue gives V* = 11.68 ± 0.11 for NGC 4874. In order to 
transform V t into M™ T , some estimation of the Coma cluster distance is required. Kelson 
et al. (2000) reported the Cepheid distance to Virgo (m — M ) = 31.03 ± 0.16. Jensen et 
al. (2001) measured the SBF distance of the Virgo galaxies NGC 4472 and NGC 4406, 
obtaining (m — M ) = 31.06 ± 0.10 and (m — M ) = 31.17 ± 0.14 respectively. We adopt the 
mean of the three previous values as our Virgo distance, i. e. (m — M)v irgo = 31.08 ± 0.14. 
Kavelaars et al. (2000) using the GCLF obtained a relative Coma-Virgo distance modulus 
of 4.06 ± 0.11. On the other hand, D’Onofrio et al. (1997), reported a fundamental plane 
relative Coma-Virgo distance of 3.55 ± 0.15. Adopting the mean of both, 3.81 ± 0.16, the 
Coma distance results (m — M) c om a = 34.89 ± 0.20. 
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The Coma cluster is located close to the Galactic Pole, with nearly negligible foreground 
extinction. Using the expresions in Postman & Lauer (1995) for the extintion and the K- 
correction, and asuming E(B — V) — 0.01 (Baum et al. 1997) and a redshift z = 0.024, 
we obtain Ar = 0.024 and Kr = 0.025 mag. After the correction of these effects, we obtain 
My° l = —23.26 ± 0.11 and hence a specific frequency Sn = 9.0 ± 2.8. 

For the remaining galaxies, SBF analysis results are presented in Table 7, in the same 
way as for NGC 4874, and for those galaxies that have high enough angular size, the radial 
structure of those GCSs are shown in Figure 8. The extrapolation to the central regions of 
all the galaxies is also shown in the first row of each galaxy in Table 7. In Table 9, My OT is 
listed in column 2. For the galaxies IC 3959, MCG +5 —31 —063, NGC 4673, and IC 3651, 
V t does not appear in the RC3 calatogue, and My OT was estimated using B t and assuming 
a mean color for these galaxies of (B — V) = 1.1. In column 3, the distance of the galaxy 
to the central galaxy NGC 4874 is given, and final results for total populations (Nqq) are 
listed in column 4 and the derived Sn in column 5. 

The distance uncertainty produces an additional error of 18% in the results of Sn- It 
must be considered when quoting Sn for a single galaxy. However, it has the same effect in 
Sn for all the Coma galaxies and does not introduce further internal dispersion in any relative 
distribution of Sn values of Coma galaxies. In Table 9, the uncertainty in the distance has 
not been included in the Sn results. 


6. Discussion and Conclusions 

In the following, we discuss first the scenarios for elliptical-galaxy formation in the light 
of the obtained results, and second, the possibility that a merging process is at work in 
Coma. 


6.1. Elliptical Galaxy Formation Scenarios 

We have measured Sn in 17 elliptical galaxies located in Coma. In Table 9 we present 
a summary of the obtained results. For each galaxy we show the value of My 0r (obtained 
from the RC3 value and adopting a distance modulus for Coma of 34.89 ±0.20), the distance 
to the central galaxy, NGC 4874, the obtained total population of GCs, and the result for 
Sn- The uncertainty in the distance has not been included in the results, it produces an 
additional error of 18% in Sn- These results reveal enormous differences in Sn among similar 
galaxies. In particular, Sn varies by an order of magnitude from galaxy to galaxy. Extreme 
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cases are: a) at the lower end, NGC 4673 has Sn = 1.0±0.4, this value being typical of spiral 
or irregular galaxies, but surprising for an elliptical galaxy; b) on the other hand, MCG ±5 
—31 —063 has Sn = 13.0 ± 4.2 and IC4051 has Sn = 12.7 ± 3.2, which are similar to the 
values found in supergiant cD galaxies, but not in “normal” elliptical galaxies. 

Blakeslee et al. (1997) also performed an SBF analysis and reported Sn results for NGC 
4874, NGC 4889 and NGC 4839. They obtained Sn = 9.3 ± 2.0 for NGC4874; 5.7 ± 1.3 for 
NGC 4889, and 4.6 ± 1.5 for NGC 4839. These results are compatible (within the error bars) 
with ours. We obtained for these galaxies Sn = 9.0 ± 2.2, Sn = 4.0 ± 1.2 and Sn = 7.0 ± 1.9 
respectively. Perhaps there is a small difference in the case of NGC 4839. Note that results 
reported by Blakeslee et al. (1997) are “ metric — Sn ”; he., Sn calculated within a radius of 
40 kpc. If GCs and halo light follow the same radial distribution, Sn does not change with 
radius and the “metric — Sn” would be identical to the global Sn- The results for NGC 
4839 probably indicate that NGC 4839 GCS is more extended than its halo. Indeed, NGC 
4839 has the most extended GCS of our sample (see Fig. 8). 

Woodworth & Harris (2000) studied HST images of 1C 4051 and proved that a central 
location in a rich cluster environment is not required to form a high population of GCs. 
They obtained a Sn equal to 11 ± 2. In our work, IC 4051, a “normal” elliptical near the 
cluster core, has a high Sn of 12.7 ± 3.2. Furthermore, NGC 4874, the central supergiant 
cD galaxy in Coma, also exhibits a relatively high specific frequency, Sn = 9.0 ± 2.2. The 
remaining galaxies studied have Sn in the range [2, 7], the mean value being Sn = 5.1. 

Why do IC 4051 and NGC 4673, galaxies with similar absolute magnitudes, have differ¬ 
ences of a factor of twelve in SV? Perhaps this is due to differences in the environment. But 
then, why does MCG ±5 —31 —063 have Sn five times bigger than IC 4041, if both galaxies 
are located on the border of the Coma cluster core and have the same absolute magnitudes?. 
In order to study possible relations between Sn and environment, Sn versus the distance R 
to the central Coma galaxy, NGC 4874 is plotted in Fig. 9. No clear trend is found in this 
plot, which suggests that Sn does not depend significantly on the environment in Coma. On 
the other hand, Fig. 10 shows Sn versus My OT of each galaxy. No relation is found between 
Sn arid My OT . The figure is completely dominated by the dispersion of the points, which is 
greater than the error bars and must therefore be real. 

Formation scenarios for giant ellipticals tend to fall into three basic classes of models: 
a) in situ models, in which the galaxy condenses by dissipative collapse of gas clouds in one 
or more major bursts, b) mergers of gas-rich systems, probably disk-type galaxies, and c) 
accretion of smaller satellites. Various combinations of these extremes are also possible. 

In situ models predict a correlation between galaxy properties and the GCS. If we assume 
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a universal efficiency of GC formation per unit total initial galaxy mass (McLaughlin 

Mgc 

M* + M gas 

where e is the efficiency parameter, M* is the mass of the visible stellar component 
galaxy, M gas is the mass of gas in or around the galactic halo, and Mgc is the mass 
GCS, then Sn must be related with the luminosity of the host galaxy according to: 

S N ~ £(1 + 

The term accounts for the systematic increase in mass-to-light ratio with galaxy 
luminosity. McLaughlin (1999) predicts the general trend of Sn with luminosity. The only 
unknown parameter affecting Sn would be the ratio of gas to stellar mass. In this paper no 
evidence is found for a relation between galactic luminosity and Sn- 

On the other hand, if a range of possible GC formation efficiencies is allowed, the merger 
model can account for the specific frequency range observed for elliptical galaxies (Ashman 
& Zepf 1992). If the merger model is to explain the high — Sn phenomenon!, GC formation 
must be more efficient in the high — Sn galaxies. 

Forte, Martinez, & Muzzio (1982) and Muzzio (1987) have suggested that in clusters 
of galaxies a population of intergalactic GCs (IGCs) should exist. A galaxy could enlarge its 
GC population via gravitational capture of IGCs. If they are distributed according to the 
gravitational potential of the whole galaxy cluster, a relation between Sn and environment is 
predicted (White 1987). The current data set shows no indication supporting this prediction. 

The fact that no single scenario seems to account for the observed specific frequencies, 
indicates the history of each galaxy should be deduced individually by suitably combining 
the models mentioned above (Woodworth & Harris 2000). To this aim, it becomes necessary 
to extend the observational information on each galaxy. Besides the specific frequency and 
the radial distribution of the GCS, it is necessary to know the details of the subpopulations 
(when these exist) and the kinematics of the GCS. In this way, detailed HST observations of 
the extreme cases mentioned above (NGC 4673 and MCG +5 —31 —063) could be a good 
starting point because these peculiar galaxies can show characteristics found nowhere else, 
and so offer valuable information in testing the different scenarios. 


1999), 

(19) 
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6.2. Subgroups and Merging in Coma 

Gurzadyan & Mazure (2001) have recently discovered the existence of three subgroups 
of galaxies in Coma, one of them associated with the cD galaxy NGC 4874 and the other two 
with NGC 4889 and NGC 4839. They conclude that the non-stationarity of the dynamical 
processes at work in the Coma core is due to the merging of small-scale groups of galaxies. 
In this context, each subgroup formed separately and then the merger between the different 
groups took place. If this scenario is valid, is there any relation between Sn and environment 
and/or My OT inside each subgroup? 

In order to analyze this question, we restricted our figures to galaxies belonging to 
Gurzadyan & Mazure (2001) subgroup 2 and studied in this paper: NGC 4874, IC 4012, IC 
4041, IC 3976, and IC 3959. In Fig. 11 Sn versus host galaxy magnitude is plotted, while 
Fig. 12 shows Sn versus the distance to the galaxy NGC 4874, which is very close to the 
center of subgroup 2. No relation between Sn and My OT is found from Fig 11. But there is 
an apparent trend in Fig. 12: Sn is bigger in high density environments. Is this trend real 
or is it an artifact of the low number of galaxies considered?. If this result is confirmed, this 
will be a strong argument in favor of the IGCs model, so the next step in this work will be 
to enlarge the number of galaxies of the study. 


We are very grateful to W. E. Harris for reading a preliminary version of this paper 
and to J. J. Fuensalida for his helpful comments. This article is based on observations made 
with the 2.5 m Isaac Newton Telescope operated on the island of La Palma by the ING 
in the Spanish Observatorio del Roque de Los Muchachos. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. This research has made use of the Digitized Sky Survey, produced 
at the Space Telescope Science Institute at Baltimore under U.S. grant NAGW-2166. This 
research has been supported by the Instituto de Astrofisica de Canarias (grant P3/94), 
the DGESIC of the Kingdom of Spain (grant PI97-1438-C02-01), and the DGLII of the 
autonomous government of the Canary Islands (grant PI1999/008). 
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Fig. 1.— Central chip of the CCD mosaic in the target field RC-1. The nine galaxies studied 
in this chip are labeled. 
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Fig. 2.— Synthetic images generated for the internal consistency test of the SBF to measure 
GCS properties. Image (a) is noise-free. Noise increases from left to right and from top to 
bottom. See text and Table 5 for details. 
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Fig. 3.— SBF analysis of the synthetic images with different noise levels. Points represent 
the power spectrum of each image, and the solid line is the SBF fitting (eq. 14). See Table 5 
for details. 
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Fig. 4.— Luminosity function of all detected objects in the NGC 4874 image binned to 0.3 
mag. Solid line represents the BGLF obtained by Tyson (1988) and scaled to our counts. 
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Fig. 5.— Luminosity functions (binned to 0.3 mag) of all detected objects in each galaxy 
studied. Solid lines represents the BGLF obtained by Tyson (1988) and scaled to our counts. 
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Fig. 6.— SBF analysis of the eight ring-shaped regions studied in NGC 4874. Points are the 
power spectrum of the masked ring-shaped region, and solid lines are the results of the SBF 
fit. See details in Table 7. 
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Fig. 7.— Radial dependence of the GCs surface density in NGC 4874. Open dots: results 
obtained with SBF in this work; Filled dots: results obtained using the HST and traditional 
techniques (Harris et al. 2000). 





Fig. 8.— Radial structure of the GCS of the nine galaxies with high enough angular size 
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Table 1. Target fields 


Target field 

a (2000) 

5(2000) 

Galaxies studied 

RC-1 

13 h 00 m 

14 s .4 

+28°00'46"9 

NGC 4874, NGC 4889, NGC 4886, 

IC 4012, IC 4021, IC 4026, 

IC 4041, IC 4045, I C4051, 

IC 3976, IC 3959, MCG +5 -31 -063 

RC-2 

12 h 57 m 

23 s .0 

+27°43'37"9 

NGC 4839, NGC 4840, NGC 4816 

RC-8 

13 h 45 m 

26 s .0 

+27°06 , 12('0 

NGC 4673 

RC-9 

12 h 40 m 

53 s . 1 

+26°43 , 39'/9 

IC 3651 
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Table 2. Summary of observations 


Target field 

Date 

Exposure (s) 

Filter 

Seeing 

RC-1 

Apr 

25 

6 

X 

300 

R 

(Sloan) 

l'.'l 

RC-1 

Apr 

27 

7 

X 

300 

R 

(Sloan) 

171 

RC-2 

Apr 

27 

12 

X 

300 

R 

(Sloan) 

171 

RC-8 

Apr 

27 

12 

X 

300 

R 

(Sloan) 

170 

RC-9 

Apr 

27 

12 

X 

300 

R 

(Sloan) 

172 








35 


Table 3. 


Zero points and colour terms of each chip of the CCD 


Zero point Colour term 


a i=24.422 6i=-0.111 

a 2 =24.671 b 2 =-0.320 

03=24.682 b 3 =-0.178 

04=24.654 b 4 =-0.132 
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Table 4. Photometric calibration 


Galaxy rrv[ 


NGC 4874 

24.44 

± 

0.01 

NGC 4889 

24.44 

± 

0.01 

NGC 4886 

24.44 

± 

0.01 

IC 4012 

24.44 

± 

0.01 

IC 4021 

24.44 

± 

0.01 

IC 4026 

24.44 

± 

0.01 

IC 4041 

24.44 

± 

0.01 

IC 4045 

24.44 

± 

0.01 

IC 4051 

24.44 

± 

0.01 

IC 3976 

24.21 

± 

0.01 

IC 3959 

24.36 

ds 

0.02 

MCG +5 -31 -063 

24.36 

± 

0.02 

NGC 4839 

24.22 

± 

0.01 

NGC 4840 

24.22 

± 

0.01 

NGC 4816 

24.37 

± 

0.02 

NGC 4673 

24.43 

± 

0.01 

IC 3651 

24.43 

± 

0.01 
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Table 5. Results of the SBF consistency test 


Experiment 

^oise /Pixel) 2 

/v in 

iV GC 

Pi (e /pixel) 2 

Po (e /pixel) 2 


(a) 

0 

10000 

2.55 

2920 ± 40 

10500 ± 140 

(b) 

625 


628 

2890 ± 60 

10400 ± 200 

(c) 

3025 

10000 

3027 

2900 ± 100 

10400 ± 400 

(d) 

4900 

10000 

4900 

2900 ± 150 

10400 ± 500 

(e) 

10000 


9994 

3000 ± 300 

10800 ± 1100 

(f) 

15625 


15568 

3500 ± 500 

12600 ± 1800 
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Tablc 6. Results for cr| G , in units of 10 “ 4 x 


Galaxy crjj| 


NGC 4874 

29.4 ± 5.9 

NGC 4889 

26.0 ± 5.2 

NGC 4886 

26.0 ± 5.2 

IC 4012 

18.0 ± 3.6 

IC 4021 

16.0 ± 3.2 

IC 4026 

17.0 ± 3.4 

IC 4041 

20.0 ± 4.0 

IC 4045 

20.0 ± 4.0 

IC 4051 

26.0 ± 5.2 

IC 3976 

12.1 ± 2.4 

IC 3959 

15.3 ± 3.1 

MCG +5 -31 -063 

13.7 ± 2.7 

NGC 4839 

10.1 ± 2.0 

NGC 4840 

10.0 ± 2.0 

NGC 4816 

17.5 ± 3.5 

NGC 4673 

14.2 ± 2.8 

IC 3651 

13.7 ± 2.7 
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Table 7. Properties of the ring-shaped regions analyzed in each galaxy and SBF results 


7* min 

(") 

Tmax 

(") 

A 

("? 

Pa x 10 4 

( e~ \ 2 (a) 

°GC x 1q4 
( e- l 2 (a) 

Ngc 

GC/(") 2 

7V^ ion 

^ sxpix ) 

^ sxpix J 

NGC 4874 

0 

12.6 

452 

- 

- 

- 

613 ± 130 b > 

12.6 

25.9 

1410 

130 ± 13 

101 ± 14 

1.36 ± 0.29 

1910 ± 410 

25.9 

37.0 

1479 

77.0 ± 3.4 

47.6 ± 6.8 

0.64 ± 0.29 

950 ± 250 

37.0 

52.5 

3781 

71.2 ± 3.3 

41.8 ± 6.8 

0.56 ± 0.17 

2130 ± 640 

52.5 

65.6 

6408 

56.9 ± 3.2 

27.5 ± 6.7 

0.37 ± 0.17 

2370 ± 900 

65.6 

102.3 

13022 

52.1 ± 2.2 

22.7 ± 6.3 

0.31 ± 0.14 

4000 ± 1700 

102.3 

125.1 

15031 

42.3 ± 1.2 

13.8 ± 6.0 

0.19 ± 0.13 

2800 ± 1600 

125.1 

141.8 

24153 

38.2 ± 3.1 

8.8 ± 6.7 

0.12 ± 0.11 

2900 ± 2600 

141.8 

161.4 

18778 

31.8 ± 1.2 

2.4 ± 6.0 

0.03 ± 0.09 

600 ± 1700 

NGC 4889 

0 

10.4 

250 

- 

- 

- 

158 ± 42 W 

10.4 

30.5 

3161 

73.1 ± 3.3 

47.1 ± 6.2 

0.63 ± 0.17 

2000 ± 540 

30.5 

45.7 

4722 

46.4 ± 3.1 

20.4 ± 6.1 

0.27 ± 0.13 

1300 ± 610 

45.7 

71.9 

12767 

40.5 ± 3.1 

14.5 ± 6.1 

0.19 ± 0.12 

2500 ± 1500 

71.9 

104.7 

10904 

37.7 ± 3.1 

11.7 ± 6.1 

0.16 ± 0.11 

1700 ± 1200 

104.7 

121.2 

6287 

38.1 ± 3.1 

12.1 ± 6.1 

0.16 ± 0.11 

1020 ± 690 

121.2 

157.1 

12558 

31.3 ± 5.0 

5.3 ± 7.2 

0.07 ± 0.11 

900 ± 1400 

157.1 

180.3 

15044 

26.9 ± 3.0 

0 

0 

0 

180.3 

222.3 

22139 

26.2 ± 2.1 

0 

0 

0 

NGC 4886^) 

0 

3.5 

45 

- 

- 

- 

26 ± 9 W 

3.5 

24.1 

1628 

80.2 ± 6.4 

42.5 ± 10.2 

0.57 ± 0.20 

931 ± 325 

24.1 

33.4 

2727 

45.0 ± 4.1 

7.3 ± 8.9 

0.10 ± 0.13 

267 ± 354 

IC 4012 

0 

6.3 

127 

- 

- 

- 

17 ± 11 W 

6.3 

39.0 

4605 

27.8 ± 3.1 

9.8 ± 4.8 

0.13 ± 0.09 

610 ± 410 

39.0 

46.7 

3692 

18.0 ± 1.1 

0 

0 

0 

IC 4021 

0 

5.0 

85 

- 

- 

- 

6±6>‘l 

5.0 

27.1 

2229 

21.0 ± 1.1 

5.0 ± 3.4 

0.07 ± 0.07 

150 ± 160 

27.1 

46.7 

6011 

16.4 ± 1.0 

0 

0 

0 

IC 4026 

0 

6.5 

129 

- 

- 

- 

78 ± 22 W 

6.5 

16.5 

694 

62.1 ± 6.1 

45.1 ± 7.0 

0.61 ± 0.17 

421 ± 120 

16.5 

28.2 

1915 

22.8 ± 3.1 

5.8 ± 4.6 

0.08 ± 0.08 

150 ± 150 

28.2 

56.7 

16810 

17.2 ± 1.1 

0 

0 

0 

IC 4041 

0 

6.4 

249 

- 

- 

- 

74 ± 30 ^ 

6.4 

25.8 

1336 

42.2 ± 3.1 

22.2 ± 5.1 

0.30 ± 0.12 

399 ± 160 

25.8 

44.3 

11113 

21.2 ± 1.1 

0 

0 

0 

IC 4045 

0 

13.7 

675 

- 

- 

- 

146 ± 74 W 

13.7 

37.3 

3654 

36.1 ± 3.1 

16.1 ± 5.1 

0.22 ± 0.11 

790 ± 400 

37.3 

52.5 

3841 

23.0 ± 3.0 

3.0 ± 5.0 

0.04 ± 0.08 

150 ± 310 

52.5 

74.7 

13400 

19.7 ± 1.1 

0 

0 

0 


IC 4051 


o 


5.5 


126 


118 ± 23 W 




















Table 7—Continued 


T min 

Tmax 

A 

P 0 x 10 4 

°GC x 10 4 

^GC 

N resi°n 

(") 

(") 

(") 2 

f e~ A 2 (a) 

i e- \ 2 (a) 

GC/(") 2 


^ sXpix J 

sXpix J 


5.5 

34.9 

4093 

95.5 ± 4.4 

69.5 ± 6.8 

0.94 ± 0.18 

3830 ± 730 

34.9 

47.2 

4367 

43.2 ± 3.1 

17.2 ± 6.1 

0.23 ± 0.12 

1010 ± 520 

47.2 

60.6 

9216 

36.4 ± 2.1 

10.4 ± 5.6 

0.14 ± 0.10 

1290 ± 920 

60.6 

96.9 

23739 

26.1 ± 2.1 

0 

0 

0 

IC 3976 

0 

8.0 

196 

- 

- 

- 

124 ± 35 W 

8.0 

16.7 

862 

43.8 ± 6.1 

31.7 ± 6.6 

0.64 ± 0.18 

550 ± 150 

16.7 

25.9 

1363 

18.2 ± 3.0 

6.1 ± 3.8 

0.12 ± 0.09 

170 ± 120 

25.9 

36.2 

1900 

14.2 ± 3.0 

2.1 ± 3.8 

0.04 ± 0.08 

80 ± 150 

36.2 

66.6 

6055 

11.5 ± 3.0 

0 

0 

0 

IC 3959 

0 

11.5 

459 

- 

- 

- 

105 ± 46 W 

11.5 

33.2 

2878 

30.2 ± 2.1 

14.9 ± 3.7 

0.23 ± 0.10 

660 ± 280 

33.2 

44.9 

2821 

18.1 ± 1.1 

2.8 ± 3.3 

0.04 ± 0.06 

120 ± 170 

44.9 

89.8 

22707 

15.2 ± 1.0 

0 

0 

0 




MCG +5 

—31 -063 



0 

1.9 

12 

- 

- 

- 

8 ± 2 W 

1.9 

15.9 

765 

55.3 ± 7.1 

41.6 ± 7.6 

0.64 ± 0.18 

486 ± 140 

15.9 

24.9 

1317 

38.1 ± 3.1 

24.4 ± 4.1 

0.37 ± 0.12 

491 ± 160 

24.9 

38.3 

2887 

23.1 ± 2.0 

9.4 ± 3.4 

0.14 ± 0.08 

420 ± 230 

38.3 

54.0 

7195 

20.5 ± 0.7 

6.8 ± 2.8 

0.10 ± 0.07 

750 ± 570 

54.0 

133.2 

43212 

13.5 ± 0.6 

0 

0 

0 

NGC 4839 

0 

31.6 

3144 

- 

- 

- 

1355 ± 470 w 

31.6 

38.2 

1730 

32.1 ± 5.0 

22.0 ± 5.4 

0.43 ± 0.15 

750 ± 260 

38.2 

60.5 

3880 

22.3 ± 1.1 

12.2 ± 2.3 

0.24 ± 0.09 

930 ± 350 

60.5 

73.6 

6452 

22.1 ± 2.1 

11.0 ± 2.9 

0.22 ± 0.09 

1390 ± 580 

73.6 

111.9 

7756 

18.2 ± 1.1 

7.1 ± 2.3 

0.14 ± 0.07 

1080 ± 540 

111.9 

138.5 

13850 

14.9 ± 1.6 

4.8 ± 2.6 

0.09 ± 0.07 

1300 ± 970 

138.5 

197.2 

13249 

13.1 ± 1.0 

3.0 ± 2.2 

0.06 ± 0.06 

780 ± 790 

197.2 

240.1 

24147 

13.0 ± 1.0 

2.9 ± 2.2 

0.06 ± 0.06 

1400 ± 1400 

240.1 

277.4 

25663 

12.1 ± 1.0 

2.0 ± 2.2 

0.04 ± 0.05 

1000 ± 1300 

277.4 

349.6 

259328 

10.3 ± 1.0 

0 

0 

0 

NGC 4840 

0 

12.9 

545 

- 

- 

- 

76 ± 38 W 

12.9 

28.1 

1705 

17.1 ± 1.6 

7.1 ± 2.6 

0.14 ± 0.07 

240 ± 120 

28.1 

44.9 

4182 

15.5 ± 1.6 

5.5 ± 2.6 

0.11 ± 0.07 

450 ± 290 

44.9 

69.2 

8811 

13.0 ± 0.8 

3.0 ± 2.2 

0.06 ± 0.06 

520 ± 530 

69.2 

108.9 

18358 

9.5 ± 0.8 

0 

0 

0 

NGC 4816 

0 

28.0 

2464 

- 

- 

- 

1590 ± 420 W 

28.0 

39.5 

2357 

60.4 ± 6.1 

42.9 ± 7.0 

0.65 ± 0.17 

1520 ± 400 

39.5 

49.6 

3053 

38.2 ± 3.1 

20.7 ± 4.7 

0.31 ± 0.12 

950 ± 370 

49.6 

61.3 

4313 

29.1 ± 1.1 

11.6 ± 3.7 

0.17 ± 0.09 

750 ± 390 

61.3 

73.2 

6025 

23.0 ± 2.1 

5.5 ± 4.1 

0.08 ± 0.08 

500 ± 480 

73.2 

104.3 

18638 

20.5 ± 0.7 

3.0 ± 3.6 

0.04 ± 0.06 

800 ± 1100 

104.3 

145.4 

34964 

17.2 ± 0.6 

0 

0 

0 

145.4 

186.1 

44838 

17.4 ± 1.1 

0 

0 

0 

















Tabic 7—Continued 


T min 

T’max 

A 

P 0 x 10 4 

°GC x 1q4 

Ngc 

n gE° U 

(") 

(") 

(") 2 

( e~ t 2 (a) 

( e ~ V (a) 

GC/(") 2 


sxpix ) 

\ sxpix J 


186.1 

241.4 

45584 

17.9 ± 1.1 

0 

0 

0 

241.4 

374.0 

193496 

16.5 ± 1.1 

0 

0 

0 

NGC 4673 

0 

15.6 

770 

- 

- 

- 

105 ± 69 w 

15.6 

30.8 

1810 

24.4 ± 4.1 

10.2 ± 5.0 

0.14 ± 0.09 

250 ± 160 

30.8 

41.3 

2370 

18.4 ± 0.7 

4.2 ± 2.9 

0.06 ± 0.06 

135 ± 140 

41.3 

58.1 

3140 

20.3 ± 1.1 

6.1 ± 3.0 

0.08 ± 0.06 

260 ± 190 

58.1 

103.5 

13364 

14.0 ± 0.6 

0 

0 

0 

103.5 

144.4 

36976 

15.5 ± 0.6 

0 

0 

0 

IC 3651 

0 

17.6 

1005 

62.9 ± 6.1 

49.2 ± 6.7 

0.66 ± 0.17 

670 ± 170 

17.6 

39.0 

3810 

32.2 ± 2.1 

18.5 ± 3.4 

0.25 ± 0.10 

950 ± 380 

39.0 

60.8 

6073 

21.0 ± 0.6 

7.3 ± 2.8 

0.10 ± 0.04 

600 ± 240 

60.8 

87.9 

12278 

16.7 ± 0.6 

3.0 ± 2.8 

0.04 ± 0.05 

490 ± 610 


a l pix = 0.333". 

b Estimated number of GCs in the central region of the galaxy. See text for details. 

c This galaxy is located in the NGC 4889 ring-shaped region with [r m i n , r maa; ] = 
[71.9 // , 104.T // ], so ctq C of NGC 4889 region has been also subtracted to Pq. 
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Table 8. 


P 2 - results for different GCLF parameters. See text for details. 


n gc 


Galaxy 

(27.22, 1.40) 

(27.62, 1.40) 

(27.42. 1.35) 

(27.42, 1.45) 

Adopted 

NGC 4874 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

NGC 4889 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

NGC 4886 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

IC 4012 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

IC 4021 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

IC 4026 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

IC 4041 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

IC 4045 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

IC 4051 

0.088 

0.051 

0.060 

0.075 

0.067 ± 0.016 

IC 3976 

0.058 

0.033 

0.039 

0.049 

0.045 ± 0.011 

IC 3959 

0.076 

0.044 

0.052 

0.065 

0.059 ± 0.014 

MCG +5 -31 -063 

0.076 

0.044 

0.052 

0.065 

0.059 ± 0.014 

NGC 4839 

0.059 

0.034 

0.040 

0.050 

0.046 ± 0.011 

NGC 4840 

0.059 

0.034 

0.040 

0.050 

0.046 ± 0.011 

NGC 4816 

0.078 

0.045 

0.053 

0.066 

0.060 ± 0.014 

NGC 4673 

0.087 

0.050 

0.059 

0.074 

0.067 ± 0.016 

IC 3651 

0.087 

0.050 

0.059 

0.074 

0.067 ± 0.016 







Table 9. Summary of results 


Galaxy 

M™ t 

^?NGC4874 0 

*GC 

Sjv 

NGC 4874 

-23.26 

± 

0.11 

0 

18200 ± 4100 

9.0 ± 2.2 

NGC 4889 

-23.45 

± 

0.11 

8.05 

9600 ± 2600 

4.0 ± 1.2 

NGC 4886 

-21.08 

± 

0.10 

7.43 

1220 ± 480 

1.8 ± 1.8 

IC 4012 

-20.00 

± 

0.11 

10.68 

620 ± 410 

6.2 ± 4.1 

IC 4021 

-20.09 

± 

0.10 

10.88 

160 ± 160 

1.5 ± 1.5 

IC 4026 

-20.35 

± 

0.15 

12.62 

650 ± 190 

4.7 ± 1.5 

IC 4041 

-20.59 

± 

0.20 

16.38 

470 ± 160 

2.7 ± 1.1 

IC 4045 

-21.00 

± 

0.10 

19.85 

1100 ± 510 

4.4 ± 2.1 

IC 4051 

-21.74 

± 

0.15 

19.95 

6300 ± 1300 

12.7 ± 3.2 

IC 3976 

-20.24 

± 

0.15 

6.78 

920 ± 250 

7.4 ± 2.2 

IC 3959 

-20.95 

± 

0.16: a 

12.67 

880 ± 330 

3.7 ± 1.5 

MGC +5 -31 -063 

-20.55 

± 

0.14: a 

8.93 

2150 ± 650 

13.0 ± 4.2 

NGC 4839 

-22.88 

± 

0.10 

43.12 

10000 ± 2500 

7.0 ± 1.9 

NGC 4840 

-21.24 

± 

0.13 

37.43 

1290 ± 620 

4.1 ± 2.0 

NGC 4816 

-22.10 

± 

0.20 

52.60 

6100 ± 1400 

8.8 ± 2.6 

NGC 4673 

-22.14 

± 

0.18: a 

217.07 

750 ± 290 

1.0 ± 0.4 

IC 3651 

-21.68 

± 

0.14: a 

290.32 

2700 ± 780 

5.7 ± 1.8 


Estimated magnitude. See text for details. 







